Structure and dynamics of colloidal depletion gels: 
coincidence of transitions and heterogeneity 
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Transitions in structural heterogeneity of colloidal depletion gels formed through short-range 
attractive interactions are correlated with their dynamical arrest. The system is a density and 
refractive index matched suspension of 0.20 volume fraction poly(methyl methacyrlate) colloids 
with the non-adsorbing depletant polystyrene added at a size ratio of depletant to colloid of 0.043. 
As the strength of the short-range attractive interaction is increased, clusters become increasingly 
structurally heterogeneous, as characterized by number-density fluctuations, and dynamically im- 
mobilized, as characterized by the single-particle mean-squared displacement. The number of free 
colloids in the suspension also progressively declines. As an immobile cluster to gel transition is 
traversed, structural heterogeneity abruptly decreases. Simultaneously, the mean single-particle dy- 
namics saturates at a localization length on the order of the short-range attractive potential range. 
Both immobile cluster and gel regimes show dynamical heterogeneity. Non-Gaussian distributions of 
single particle displacements reveal enhanced populations of dynamical trajectories localized on two 
different length scales. Similar dependencies of number density fluctuations, free particle number 
and dynamical length scales on the order of the range of short-range attraction suggests a collective 
structural origin of dynamic heterogeneity in colloidal gels. 

PACS numbers: 82.70.Dd, 64.70.Ja, 82.70.Gg 
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I. INTRODUCTION 



Improved understanding of the origin and conse- 
quences of slow dynamics in disordered colloidal suspen- 
sions impacts their application in both established [lj and 
emerging Q technologies. Slow dynamics in colloidal gels 
and glasses can arise under a broad range of suspension 
packing fractions, potential interactions, and structures. 
At dense volume fractions, without any clear structural 
transition Q, colloidal hard spheres form a glass due to 
caging that is characterized by dynamics [4j,|5j several or- 
ders of magnitude slower than its liquid precursor. Gels 
comprised of colloids interacting through short-range at- 
tractive forces share this dynamic arrest 6]; however, 
here the mechanistic origin is unclear, since observed gel 
structures are highly variable - ranging from the rami- 
fied fractal clusters of dilute gels formed from aggrega- 
tion into deep potential wells 7] to the nearly equilibrium 
liquid structures of densely packed attractive glasses [g • 
Understanding the effect of manipulating control param- 
eters on microstructure and microdynamics is the key to 
deeper understanding of gel phase behavior, rheological 
response and processing characteristics. 

As volume fraction is increased, various phenomcno- 
logical regimes have been observed including fractal 
cluster gels |p , mobile and immobile cluster phases 
|ToL [Til fl2l Il3| . concentrated attractive gels and attrac- 
tive glasses |8|. Although open questions in all these ar- 
eas remain, the range of intermediate particle loading is 
poorly understood in particular. In this regime, typically 
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characterized by both intermediate [T^ and long-range 
|lol| structural heterogeneity, neither the micromechani- 
cal ideas based on fractal dynamics [lql nor the excluded 
volume caging mechanisms of glasses |5( can be applied. 

This regime of concentrated gelation exhibits complex 
structural transitions [HE ElES EH- In addition 
to network gels in various degrees of arrest and com- 
paction |lfj|. ex per iments have recently revealed clus- 
ter fluid phases [101 [l2j and ordered tetrahedral chains 
[22| . Moreover, within this intermediate volume fraction 
range, clustering and gelation boundaries may lie in close 
proximity to metastable liquid /cr ystal and liquid/liquid 
equilibrium phase boundaries |23| . Gel boundaries also 
can display a strong kinetic component that leads to tran- 
sient gelation in non-density matched suspensions [24|. 

Mechanistic descriptions of the dynamic transitions 
that accompany these regimes include predictions of par- 
ticle localization due to strong short-range interactions 
from direct simulation (2^ , idealized mode coupling the- 
ory , and naive mode coupling theory |27| . Simulation 
and theory have attributed the arrest that leads to dy- 
namical heterogeneity to effective particle-particle bonds 
arising from short-range attractive potentials |2^| , to the 
development of entropic barrier hopping or to the 
emergence of two distinct particle populations with long 
exchange times 30] . Earlier work has suggested a role for 
arrested spinodal decomposition and other equilibrium 
transitions in the development of ultimate gel structure 
|20| . Recent studies have noted the possible perturba- 
tive effect of metastable equilibrium phase boundaries 
|l4j . and their sensitivity to the presence of long-range 
charge repulsions that are common in colloidal suspen- 
sions jU . Each mechanistic description yields particular 
implications for gel boundaries, structure, kinetics, and 
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viscoelasticity. Thus, the relative role and interaction of 
these mechanisms in model colloidal gel systems should 
be carefully assessed through experiment. 

Because of their availability, most quantitative com- 
parisons to date between theory and simulation on the 
one hand and experiment on the other have focused on 
macroscopically measurable quantities such as gel bound- 
aries and collective behavior such as the dynamic struc- 
ture factor. Previously, pair dynamics of particles in 
gel networks have been quantified |32T |. To further dis- 
criminate between the competing mechanisms suggested 
by very recent theory and simulations [3(j, comparison 
to microscopic measurements of localization length, dy- 
namical heterogeneity, structural populations, and bar- 
rier hopping are required. Suggestion that rare dynam- 
ical processes play a role in gel response (2^, |33 adds 
another argument for measurement of single particle dis- 
tributions of dynamics. It is the aim of this work to fill 
these gaps by direct quantification of local gel structure 
and dynamics through confocal microscopy. In particu- 
lar, Lu et al and Sedgwick et al |l2j have recently 
reported complex structural transitions in this regime. 
The next step in this area should be to correlate struc- 
ture to dynamics. 

Our gel system takes advantage of the well- 
characterized nature of the polymer depletion inter- 
action. We choose monodisperse sterically stabilized 
poly(methyl methacrylate) spheres |3J| that can be dis- 
persed in density and refractive index matched solvent 
mixtures to minimize the effects of polydispersity, sedi- 
mentation, and van der Waals forces, which complicate 
experimental observation and modeling of gel behavior. 
Depletion is chosen from a range of possibilities |33l l35| 
to provide an attractive force because it has easily tun- 
able experimental parameters controlling both range and 
strength of attraction and it has been well described by 
theory [M El III] . We vary the strength of the deple- 
tion potential to determine the structural and dynamic 
dependence for a given attractive range. It has recently 
been reported that the range of attraction affects the 
cluster morphology [lOf, Thus we fix this parameter at 
about four percent of the radius of the colloid in a regime 
where stable, tenuous gels are expected and fine struc- 
tural and dynamic changes can be differentiated. Flu- 
orescent particles in refractive index matched solvents 
lend themselves to direct observation via confocal mi- 
croscopy |32| yJI , which allows local material properties 
to be probed. 

This paper focuses on four measures that quantify 
both structural and dynamical heterogeneity in gels. For 
structure, we report short-range structure by contact or 
bond number distributions and long-range structure via 
number density fluctuations. For dynamics we measure 
single particle mean-squared displacement as well as the 
van Hove self-correlation function. Here, we focus on 
the effects of depletion pair potential strength on the 
quantities. We find rich variation in, and correlation of, 
structural and dynamic measures as the strength of pair 



interactions is increased. 



II. METHODS AND MATERIALS 

A. Synthesis and characterization of monodisperse, 
density matched colloidal system 

The method to synthesize the poly(12-hydroxystearic 
acid) (PHSA) stabilized monodisperse poly(methyl 
methacr ylat e) (PMMA) colloids is as reported elsewhere 
[34l l39l liol ] . The fluorescent dye incorporated in the 
dispersion polymerization of the methyl methacrylate 
monomer is Nile Red (X ex = 548 nm, \ em — 567 nm) [34| . 
The PHSA stablizer is covalently bound to the PMMA 
colloids in a locking step. After synthesis, the particles 
are washed 8 times in hexane, dried for at least 48 hours, 
characterized, and dispersed in solvents for experiments 
by vortex mixing and sonication. 

By scanning electron microscopy of the dry particles, 
the PMMA colloid diameter is 2.04 ^m ± 3.83% (Fig. 
H^a)). By differential centrifugation, we determined a 
density matching solvent composition to be 65.6% cy- 
clohexyl bromide (CHB, 98% purity, Sigma-Aldrich) and 
34.4% decalin (99+% mixture of cis and trans, Sigma- 
Aldrich). At these conditions, solvent equilibrated par- 
ticles do not sediment after ten minutes of centrifuging 
at 1500 rpm (240 g). Thus, ^ < 6 x 10" 3 . The par- 
ticles and solvent mixture are also nearly refractive in- 
dex matched at these conditions. The dried particles 
are stably redispersed in this solvent composition. The 
narrow polydispersity of the colloids is demonstrated by 
their quiescent crystallization, as reported at <fi = 0.49 
in Fig. ^b). The volume fraction range for fluid-crystal 
coexistence in the system is shifted relative to that for 
hard spheres. This shift is consistent with the presence 
of charge effects |4l l El li^ . 



B. Charge interaction measures 

The electrophoretic mobility was measured in a dilute 
(~1 vol %) sample of PMMA, suspended in a density- 
matched mixture of cyclohexyl bromide and decalin at 
25° C using Phase Analysis Light Scattering (Malvern 
Instruments Zetasizer Nano, ZEN 3600). Thirty experi- 
mental runs were performed on the sample. The resulting 
value of the ^-potential was obtained by averaging the 
mean values of all runs. At these solvent conditions we 
measured the viscosity of our solvent mixture as 2.06 cp 
(Schott-Gerate AVS 350 viscometer) and used the value 
of the material dependent dielectric constant, e r = 5.6, 
as reported in literature for similar ratios of the same sol- 
vents 0|- The colloid ^-potential was found to be +27 
mV. 

Because our system had no salt added, the ion con- 
centration was assumed to come from the dissociation 
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FIG. 1: Characterization of poly(12-hydroxysteric acid) sta- 
bilized poly(methyl methacrylate) colloids by (a) scanning 
electron micrographs of dry particles and (b) confocal mi- 
crographs of system crystallization at <j> = 0.49 dispersed in 
the solvent mixture CHB/decalin. 



of the ions within the organic solvent mixture. To de- 
termine the conductivity of our solvent mixture (73.9% 
mass fraction CHB), we used the estimated conductivity 
reported previously for the no-salt case of our system [45j . 
At 25°C, we obtained a conductivity of 125.7 pS/cm, by 
Walden's rule, implying an ion concentration of 3.0 x 10~ 9 
mol/L We calculate the Bjerrum length, Ah = 10.2 
nm; the Debye double layer thickness, = 1.46 /_tm; 
and the surface charge, Q e = 165 e/colloid. 



C. Depletion pair potential interaction due to 
non-adsorbing polymer 

The non-adsoring polymer used to generate the 
depletion pair potential interaction is monodisperse 
polystyrene (Pressure Chemical, Pittsburgh, PA) of 
molecular weight M w = 9.0 x 10 5 g/mole and polydis- 
persity M w /M n < 1.10. At the theta condition, the ra- 
dius of gyration (R g ) of this polymer would be 22 nm 
|46j . However, we find that the density matching mix- 
ture of CHB/decalin is a good solvent for polystyrene. 
This characterization was made by two independent mea- 
surements of solvent quality: mixed solvent static light 
scattering and intrinsic viscometry. Characterization of 
R g and overlap concentration, c* of the non-adsorbing 
polymer is important because these quantities determine 
the range (R g ) and magnitude (c*) of the depletion in- 
teraction. 

The R g of polystyrene in the decalin/cyclohexyl bro- 
mide mixture was characterized by static light scatter- 
ing (DAWN EOS, Wyatt Technology Corporation, Santa 
Barbara, CA) and application of a Zimm plot and the 
mixed solvent theory of Yamakawa ^3] • We determined 
R g = 41 ±4 nm. Applying the relationship c* = 4 ^i^ A , 

we compute the overlap concentration c*sls = 0.0054 
g/mL. 

The intrinsic viscosity, [77], is an indirect measurement 
of coil expansion and solvent quality. Futhermore, c*[ry] 
~ 1 01- By means of capillary viscometry of 6 dif- 
ferent polymer concentrations (Schott-Gerate AVS 350 
viscometer with Ubbelohde capillaries No. 52503 and 
52501) we find [77] = 190 ± 0.4 mL/g at T = 22°C. 



From this value, c* 
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0.0053 g/mL. The good agree- 



ment between these two independent measurements of 
solvent quality and coil expansion indicates that the den- 
sity matched mixture of CHB/decalin is a good solvent 
for polystyrene. The difference between our measured 
value and the theta estimate is likely related to the need 
to explicitly account for the mixed solvent quality effects 
on coil expansion as per Yamakawa |47|. 

In summary, based on polymer and colloid characteri- 
zation, we find £ = R g /a w 0.043 and c* « 0.0053 g/mL. 
For reference, the depletion potential strength used for 
experiments is listed in Table [I] in terms of the non- 
adsorbing polymer volume fraction, c/c*. For conve- 
nience, we also report the contact pair potential of all 
experimental conditions, estimated u sing the well-known 
Asakura Oosawa model of depletion [361 13?| . More pre- 
cise determination of the pair potential would be possible 
by application of, for example, PRISM [49j . 



D. Sample preparation 

Dry particles were dispersed in the density matched 
CHB/decalin solvent mixture by ~ 30 s vortex mixing 
and then ~ 10 minutes sonication. The dispersion was 
quiescently equilibrated for 10 - 12 hours. During this 
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TABLE I: Conditions of dimensionless polymer concentration and corresponding attractive potential between two particles 
in contact as calculated by the Asakura - Oosawa (AO) model of depletion for all experiments. Braces denote calculations 
formally outside of the range of the model. To calculate the corresponding c/c* or ILAo/kT for this system with a theta solvent, 
multiply by the value 0.15. 



c/c* 


0.15 


0.21 


0.26 


0.31 


0.33 


0.37 


0.41 


0.46 


0.49 


0.64 


1.03 


1.54 


UaoAT 


7.1 


9.5 


11.8 


14.2 


15.4 


17.1 


19.0 


21.3 


22.5 


29.6 


{47.4} 


{71.1} 



time, the colloid size evolves due to its interaction with 
the solvent. From measurement of contact conditions of 
the radial distribution function, g(r), in fully gelled sam- 
ples, the ultimate colloid diameter in the CHB/decalin 
mixture is estimated to be 1.9 /im. 

To produce suspensions with short-range attractive in- 
teractions leading to gelation, the particle dispersion was 
mixed for ~5 s on a vortex mixer with a stock solution 
of non- adsorbing polymer to achieve the c/c* values of 
Table 0] Specimens for CLSM visualization were formed 
by dispensing 180 fiL total volume in capillaries with an 
inner diameter of 6 mm. This volume corresponds to a 
height to diameter ratio of 1.06. The sample was covered 
by approximately 250 /iL of dcionizcd water, sealing the 
capillary to prevent evaporation of the solvent. 

After preparation, the sample was left undisturbed on 
the microscope stage for a minimum of 200 minutes prior 
to imaging. The kinetics of the gelation process varied 
between the samples, with the behavior being primar- 
ily dependent on the concentration of depletant. Gener- 
ally, in those samples where an immobilized gel structure 
was formed, the space-spanning network formed during 
the first 60 minutes after the dispersion was prepared. 
Various structural rearrangements of the gel structure 
occurred during the following 100 minutes, during which 
time some of the particles and clusters within the network 
moved. Once this rearrangement had occurred, the sam- 
ples stabilized to a steady structural state. We waited an 
additional 40 minutes before executing experiments. 



E. Confocal microscopy and image processing 

Gel specimens were visualized by confocal laser scan- 
ning microscopy (Leica TCS SP2) on an inverted micro- 
scope (Leica DMIRE2) using a 1.4 numerical aperture 
lOOx oil immersion objective and an excitation wave- 
length of 543 nm. At these conditions the optical resolu- 
tion is ± 240 nm in the plane of the objective and ± 550 
nm in the axial direction perpendicular to the objective 
plane. Here the objective plane is parallel to the cover 
slip that forms the lower surface of the capillary. Struc- 
ture and dynamics were quantified for all gel specimens. 
For structure measurements, image volumes consisting 
of 247 slices at 512x512 resolution were acquired. The 
pixel resolution was 146.4 nm in the x and y directions, 
and 122.1 nm in the z direction. Thus, the total image 
volume size was: 75x75x30 /ito 3 . To avoid wall effects 
the bottom edge of the volume was located at least 30 



/im from the cover slip. For studies of particle dynamics, 
time series of 500 two-dimensional 512x512 images (75 
x 75 fj,m 2 ) were collected. The time step between im- 
age acquisition was 0.848 s for a total duration of 424 s. 
Location and tracking of particles in the image volume 
and time series is accomplished by means of algorithms 
described by Crocker and Grier |5(j • The steps include a 
Gaussian convolution of the image pixels within a local 
region approximately the size of the particle radius, se- 
lection of the locally brightest pixel as a particle centroid, 
subpixel adjustment based on the intensity moments of 
the local region, and reconstruction of either a three di- 
mensional volume for structural measures or trajectories 
in two dimensional time series. 

The fidelity of the particle location algorithms of the 
image processing was assessed by the good correspon- 
dence between colloid volume fractions computed from 
the added mass of colloid and from the number of parti- 
cles located within an image volume for image volumes 
obtained in the regimes with high overlap concentrations 
of polymer and relatively homogenous density (average 
relative error = ± 4.1%). Furthermore, identified particle 
centroids were statistically uncorrelated with the location 
of the pixel grid. By applying the methods of Savin and 
Doyle |5l|, we determined the static error in particle lo- 
cation to be ± 35 nm (± 0.24 pixel) in the objective plane 
and ± 45 nm (± 0.37 slice depth) in the direction per- 
pendicular to the objective plane. Finally, the fidelity of 
the particle trajectories extracted from the image time 
series was assessed by comparing the self-diffusion coef- 
ficient of a 1 vol% colloidal suspension extracted from 
measurement of the colloid mean-squared displacement 
to the predicted self-diffusivity of a free particle given 
the measured solvent viscosity. This experiment, con- 
ducted in dioctyl phthalatc, yielded a relative error in 
the self-diffusion coefficient of ~ 0.1%. 



F. Quantitative measures 

We quantify sample structure by means of contact 
number distributions and number density fluctuations. 
Sample dynamics are measured by means of the mean- 
squared displacement and van Hove self-correlation func- 
tion. These measures are briefly illustrated for a strong 
gel formed at c/c* = 0.64. A sample CLSM image of this 
specimen is shown in Fig. [21 a). 
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FIG. 2: Quantitative methods to characterize system structure and dynamics shown for c/c* =0.64: (a) 2D confocal image, 
(b) contact number probability, (c) normalized number density fluctuations, and (d) van Hove self-correlation function. 



1. Contact number distributions 



2. Number density fluctuations 



The effect of depletion-induced gelation upon the local 
structure is determined by analysis of contact or bond 
number distributions. Two particles are considered to 
be in contact when the distance between their centers is 
less than the distance to the first minimum in the pair- 
correlation function, g(r) [2^J . We use this method in 
an attractive heterogeneous environment to characterize 
local variation. The first observed minimum in g(r) of the 
c/c* = 0.64 sample, which is considered to be a strong 
gel, at r/2a = 2.3 is applied to all samples for determining 
contact number. 



To find the mean contact number and variation for 
each polymer concentration, this measurement was per- 
formed on six image volumes and the average of the six 
mean contact numbers is reported. Figure^b) shows the 
distribution of contact number within one sample image 
volume. 



Within a sample volume, the long-range structure of 
a colloidal dispersion can be quantified by the mean- 
squared fluctuation of the number of particles in the sys- 
Since this quantity is proportional 



tern, 



<n-*>-<n> 2 



<N> 



to the isothermal compressibility of the sample, is re- 
lated to the integral over all space of the pair-correlation 
function g(r) in the macroscopic limit, and can be re- 
lated to the static structure factor in the low angle limit, 
S(0), it can be used to show the change in the long-range 
struct ure|l5|. A maximum in the measure of number 
density fluctuations is indicative of a preferred cluster 
size. The mean square fluctuation of the number of par- 
ticles in the sample was calculated by dividing the six in- 
dependent 75x75x30 /im 3 image volumes into cubic boxes 
of linear dimension L. L was systematically varied from L 
~ a to L ~ the z-dimension of the image volume. Figure 
IHc) shows sample number density fluctuation measures 
for c/c* = 0.64. 
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3. van Hove self-correlation function and mean-squared 
displacement 

The self part of the van Hove correlation function, 
Gs(x, t), measures the single particle displacement prob- 
ability distribution 0, |E2 ■ In one dimension: 



G s (x,t) = - 



*5[x + xi(t = 0) -Xi(t)] 



(1) 



where N is the number of particles and Xi(t) is particle i's 
position at time t. A free particle experiencin g B rownian 
motion has a Gaussian form for Gs{x,t) |52Ll53|: 



G 9 s (x,t) 



2ir{x 2 {t)) 



exp 



2<z»(t)> 



(2) 



where (x 2 (t)) is the one dimensional meansquared dis- 
placement. Dynamic measurements reported in this pa- 
per are the average of 4 separate time series in different 
locations in the sample and error bars represent the stan- 
dard error of the mean. Figure Hfd) shows Gs(x,t) at 
fixed c/c* for several instances of time. Error bars are 
plotted for every fifth data point for clarity on the log 
scale graph. 

The one dimensional mean square displacement is the 
second moment of Gs(x,t) or ((x(t) — a;(0)) 2 ) averaged 
over all particles in the system and is a collective measure 
of dynamics. For a free particle, this value grows linearly 
with time as per the Stokes-Einstein equation. Alterna- 
tively if a particle is locally bound to a gel network, its 
localization length can be derived from the plateau in 

<* 2 (i)> m. 



III. RESULTS 

Fig. |3| reports representative 2D CLSM images of 
structure as the magnitude of short-range interactions are 
increased. At the lowest concentration shown (Fig. OJa) , 
c/c* = 0.15), a mobile, dispersed liquid is observed. As 
attraction is increased, mobile clusters form (Fig. %Rb), 
c/c* = 0.21). The clusters coexist with free particles. 
As polymer depletant concentration is further increased, 
clusters become increasingly immobilized (Fig. GJc), c/c* 
= 0.26) until a bicontinuous network of clusters and voids 
is formed (Fig. E[d)-(g), c/c* = 0.31 -0.41). At these lev- 
els of pair colloid attraction, large-scale motions in the 
suspensions are no longer apparent. At c/c* ~ 0.46, a re- 
gion of maximum cluster/void heterogeneity is achieved. 
As depletant concentration is further increased (Fig. Ep), 
c/c* = 0.49), structures increasingly appear more homo- 
geneous and branched. At the maximum concentrations 
studied (Fig. EOj)-(l), c/c* = 0.64 -1.53) chain like struc- 
tures are apparent locally; however, globally the struc- 
ture is quite homogeneous. 

To quantify the local connectivity of the system, we 
compute contact number distributions for all c/c* stud- 
ied as per the methods of Section Til F II A color-coded 



representation (generated from RasMol, U-Mass) of thin 
(~ 5 /im) projections from the sampled image volumes 
is shown in Figure 0] for four different pair attraction 
strengths. Fig. 01a) shows the system at c/c* corre- 
sponding to the immobilized cluster regime. The cluster 
structure is apparent from the local regions of high con- 
tact number (orange-red coloring) . Mobile particles, with 
low contact probability, surround the clusters (purple- 
blue coloring) . The spatial extent of high contact number 
regions is highly variable. In Fig. EJb) and (c) we see that 
local contact number increases in the mean. That is, the 
number of low-contact number colloids progressively de- 
creases with most being located on the surfaces of voids. 
At the highest c/c*, shown in Fig. Efd), less spatial vari- 
ability in the contact number is apparent. Furthermore, 
few regions of very high contact number are found. 

These trends are apparent in the full contact number 
distributions plotted in Fig. ISJa) for the range of deple- 
tant concentration studied. Data are offset for clarity. 
Because of the offset, we also plot the probability of free 
particles (contact number = 0) for all the samples. These 
data represent the points on the ordinate of Fig. |5fa). 
Finally, the mean contact number of each pair attraction 
strength is plotted in Fig. ISJc). The mean contact num- 
ber passes through a maximum at c/c* = 0.37, where 
each particle is surrounded, on average, by 4.6 other par- 
ticles. The three figures quantify three qualitative ob- 
servations: First, low and high pair attraction strength 
results in nearly the same mean contact number. Second, 
intermediate pair attraction strengths (c/c* ~ 0.4) yield 
broad contact number distributions with relatively high 
mean values. Third, the number of free particles with 
zero contact number progressively decreases. As many 
as 10% of the colloids coexist outside aggregated clusters 
at c/c* ~ 0.2. By c/c* ~ 0.4, the number of zero-contact 
number colloids detected has decreased by more than an 
order of magnitude. 

The Fig. [31 images show that the amount of long- 
range heterogeneity in the sample varies dramatically as 
the strength of pair attractive interactions is increased. 
By varying the control parameter a/L, we assess the 
length scale dependence of number density fluctuations, 
per the method of section III F 21 and therefore the de- 
gree of heterogeneity can be quantified. Because of the 
non-monotonic variation of number density fluctuations 
with c/c*, data are plotted on two panels, Fig. EIa) and 
(b). At low c/c*, number density fluctuations are rel- 
atively small and insensitive to the inverse length scale 
a/L. For comparison, the compressibility of a hard sphere 

0.21 in a macroscopic limit 



suspension is 



(N) 



|54j . For c/c* > 0.35, heterogeneity rapidly increases. 
The increase in compressibility is reminiscent of spinodal 
decomposition: if there is a characteristic length scale in 
the system, it has diverged beyond ~ 30 colloid radii. 
In this intermediate concentration range, number den- 
sity fluctuations are more than two orders of magnitude 
greater than expected for a hard sphere fluid. 

The maximum in number density fluctuations is en- 
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FIG. 3: Representative 2D confocal micrographs of system structure as a function of increasing depletion pair potential strength, 
(a) c/c* = 0.15 (b) c/c* = 0.21 (c) c/c* = 0.26 (d) c/c* = 0.31 (e) c/c* = 0.33 (f) c/c* = 0.37 (g) c/c* = 0.41 (h) c/c* = 
0.46 (i) c/c* = 0.49 (j) c/c* = 0.64 (k) c/c* = 1.03 (1) c/c* = 1.54. 



countered at c/c* = 0.46. As the strength of pair in- 
teractions is increased further, Fig. Etb) reveals two in- 
teresting features of gel structure. First, the magnitude 
of number density fluctuations collapses to values of ap- 
proximately unity. Second, a characteristic length scale, 
revealed as a maximum in the measure emerges at a/L 
~ 0.1. This characteristic length scale is a clear indica- 
tion of the emergence of a gel structure with a heteroge- 
neous, string-like structure on intermediate length scales 
of about 10 colloid radii, but with a relatively homoge- 
neous structure on longer length scales. The full sequence 
of structural states is captured by Fig. Etc), which shows 
how heterogeneity evolves in the system with increasing 
pair potential strength at the characteristic inverse length 
scale a/L ~ 0.1. 

We now assess the implications of the structures re- 
ported in Figures 3-6 on system dynamics. Using the 
method of Sec. Ill F 31 we first show the effect of c/c* 
on the van-Hove self-correlation function in Fig. [7| Here 



data are plotted for one characteristic time (t = 8.48 s). 
Data plotted for other times display the same trends. 
Figure quantifies three trends in the distribution of 
single-particle displacements in the system: First, at 
low c/c*, displacement distributions are nearly consis- 
tent with the Gaussian function, consistent with the ob- 
served liquid- like structure (Fig. Eta)). Second, as c/c* is 
increased, single particle dynamics become increasingly 
localized. The trend in c/c* is monotonic. The great- 
est probability for displacement is much less than a col- 
loid radius. Third, as the strength of pair-attractions 
increased, the distribution of displacements appears in- 
creasingly non-Gaussian. 

Figure El shows information complementary to Fig. \7\ 
here at two characteristic c/c* we explore the time de- 
pendence of the distribution of single-particle displace- 
ments. Fig. Et a ) reports data at low pair attraction 
strength (c/c* = 0.26), characteristic of an immobilized 
cluster structure. Figure Etb) is for a gel (c/c* = 1.54) 
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FIG. 4: (Color online) Color coded representation depicts spatial dependence of contact number on depletant concentration for 
(a) c/c* = 0.26, (b) c/c* = 0.37, (c) c/c* = 0.41, (d) c/c* = 1.03 where warmer red/orange colors correspond to high contact 
numbers and cooler blue/purple colors correspond to low contact numbers. Free particles (contact number equals zero) are 
colored gray. 



with string-like heterogeneity. Both samples show non- 
Gaussian dynamics. At c/c* = 0.26, the Gaussian distri- 
bution appears to underpredict the distribution at both 
very short and very long displacements. It overpredicts 
the measurements at intermediate times. Although there 
is some time-dependence in the single-particle displace- 
ment distributions at c/c* = 0.26, in the strong, string- 
like gel regime, the dynamics are nearly time indepen- 
dent. 

Figures and [S] indicate the following interesting be- 
havior that warrants further analysis: First, particles are 
increasingly localized as c/c* is increased. This effect 
can be assessed by investigating the time-dependence 
of the single-particle mean-squared displacement (Fig. 



EJa)). Second, displacement distributions are highly non- 
Gaussian. This effect can be assessed by investigating 
the deviations of Gs(x,t) from a Gaussian function with 
variance equal to the measured mean-squared displace- 
ment (Fig. irUT a)). We conclude this section by discussing 
these two dynamical observations. 

In Fig. Ufa), the time dependence of the one- 
dimensional mean-squared displacement, (x 2 (t)) is plot- 
ted for all the c/c* conditions. Except for the sample 
with mobile cluster structure (c/c* = 0.19), the behav- 
ior at all c/c* is sub-diffusive. Mobile clusters (c/c* = 
0.19) approach the long-time self diffusivity of the col- 
loid volume fraction <!> = 0.2 when Stokes-Einstein dif- 
fusion coefficient is modified to include deviation from 
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FIG. 5: (a) Contact number distributions for each depletant concentration. Data are offset for clarity, (b) Probability of free 
particles (contact number = 0) for each sample, (c) Mean contact number distribution for each c/c*. 




a/L a/L ""'«*'' 

c/c 

FIG. 6: Number density fluctuation measures show the normalized mean-squared fluctuations of particle number as a function 
of sub-system size considered as it varies with attraction strength for (a) c/c* < 0.46 and (b) c/c* > 0.49. (c) Mean number 
density fluctuations for all attraction strengths for inverse length scale a/L = 0.101. 



the dilute limit caused by hydrodynamic interactions as 
per [55| . Remaining deviation is due to the presence of 
clusters. Within the temporal range probed by the ex- 
periment, the average particle mobility monotonically de- 
creases with increasing c/c*. The behavior is also increas- 
ingly sub-diffusive, until, under conditions that corre- 
spond to strong gels with string-like structure, the mean- 
squared displacement is nearly independent of time. 

Interestingly, for the strong gel structures, the local- 
ization in (x 2 (t)) is roughly comparable to the square of 
the polymer radius of gyration, R 9 , in good accord with 
the naive/PRISM mode coupling theory prediction |27| . 
The significance of the polymer radius of gyration here 
is that it sets the range of the short-range attractive pair 
potential between the colloids. To fully explore this com- 
parison further, in Fig. Efb) we plot the plateau value of 
(x 2 (t)) (taken at t = 42.4 s) scaled on the square of the 
polymer radius of gyration. In the immobilized cluster 
regime (c/c* < 0.26), the localization length decreases 



rapidly. For stronger gels {x 2 (t)) /K 2 ~ 5, which corre- 
sponds to a localization length of 2.24R g . For compari- 
son, Ref. [22I predicts r ioca/lzatlon /R g ~ 2 for a compara- 
ble system. This first direct experimental measurement 
of localization length is significant because it is predicted 
to be a principle determinant of the linear viscoelastic 
response of gels [2^. 

By empirically correlating the low c/c* and high c/c* 
localization behavior, as shown by the curves in Fig. 
EJb), we can estimate the point of the dynamic transi- 
tion to highly localized behavior. We find the transition 
c/c* = 0.42. It is instructive to compare this result to our 
earlier assessment of structural transitions based on the 
contact number (c/c* = 0.37) and number density fluctu- 
ations (c/c* = 0.46). The agreement of these structural 
and dynamic transition measures suggests a key correla- 
tion between structure and dynamics in the system. 

To estimate deviations of single-particle dynamics from 
Gaussian behavior, we adopt the method of Kob et al 
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Ax (\im) 

FIG. 7: van Hove self-correlation function for all polymer con- 
centrations at t = 8.48 s compared to the Gaussian estimate 
for a free particle. 

[5^| developed to analyze glassy dynamics. The method, 
shown in Fig. HOf a) for a gel at c/c* = 0.33, is a quan- 
titative way to assess dynamical heterogeneity in sys- 
tems with slow dynamics. This method clearly shows 
three regimes of behavior, as was qualitatively apparent 
in the original data of Figs. 0and|Hl At small and large 
displacements, positive deviations of the data from the 
Gaussian form are observed. At intermediate displace- 
ment, negative deviations from the Gaussian form are 
observed. Thus, we can identify two length scales that 
demark the regime transitions, xi* and X2*, as defined in 
Fig. HOf a). The hrst, xi*, is qualitatively similar to the 
localization length discussed in Fig. EJb). The second, 
X2*, is a new scale suggested by the data. It identifies 
a small population of highly mobile particles, consistent 
with the existence of dynamical heterogeneity in gels. 

IV. DISCUSSION 



t = 0.848 s 




Ax (|im) 



FIG. 8: van Hove self-correlation function at (a) low [c/c* 
= 0.26] and (b) high [c/c* — 1.54] depletant concentrations 
plotted for a range of times. 



The static structure of gels formed from strong short- 
range interactions at low and intermediate volume frac- 
tion has recently generated intense interest. In addition 
to the strength of the pair potential, structures observed 
are a sensitive function of potential range [lfj, charge 
[l2|. and density matching |l2fl. Our observations of a 
sequence of stable fluid, mobile cluster, immobilized clus- 
ter, and gel phases as pair potential strength is increased 
are broadly consistent with these recent reports. Fixed 
parameters in our experiments are that the colloids are 
charged (C = 165 e/colloid, n~ x = 1.46 /mi), the attrac- 
tive potential is short ranged (£ = 0.043), and the system 
is density-matched (Ap/p ~ 10~ 3 ). Our quantification of 



long-range number density fluctuations, short-range con- 
tact number, and single-particle dynamics allows us to 
buil d up on the phase diagrams available in Refs. 01 
and 

The difference between development of stable mobile 
and immobile cluster phases has been attributed to the 
degree of density matching . While our results do not 
directly test this hypothesis, the combination of structure 
and dynamical characterization reported here identifies 
the following additional features of these phases. 

First, Fig. Ela) demonstrates that mobile and immo- 
bile cluster phases can be distinguished by the behavior 
of the single-particle mean-squared displacement: mo- 
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: • c/c* = 0.21 » c/c* = 0.37 + c/c* = 0.64 
fc. c/c* = 0.41 

-♦ c/c* = 0.31 J c/c* = 0.46 eb c/c* = 1.54 




-|0 I 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

C/c* 

FIG. 9: (a) (Color online) Mean-squared displacement as a 
function of increasing c/c* relative to the theoretical mean- 
squared displacement of a particle in the equivalent medium 
with no non-adsorbing polymer, (b) Localization length es- 
timated from mean-squared displacement plateau decreases 
with increasing attraction point until it stabilizes at about 
2.2 Rg and is no longer sensitive to concentration. 



bile systems (such as c/c* = 0.21) show diffusive behav- 
ior while immobile phases (0.26 < c/c*< 1.54) display a 
plateau in the mean-squared displacement. This quan- 
titative metric complements previous observational dis- 
criminations and demonstrates that both types of phases 
can be observed at fixed density matching and po- 
tential range . 

Second, unlike the gel phase, where system dynamics 
are nearly independent of pair attraction strength (c.f. 



Fig. E2b)), the immobile cluster phases dynamics are ex- 
tremely sensitive to pair attractions. For example, as the 
immobile cluster phase envelope is traversed, the plateau 
mean-squared displacement changes by about an order 
of magnitude. 

Third, the sensitivity of dynamics to pair attractions in 
the immobile cluster phase is correlated with significant 
changes in structure. Locally, the number of free parti- 
cles changes from about 10% of the sample to about 1% 
of the sample within the phase boundaries. In addition, 
long-range density-fluctuations rapidly increase with pair 
attraction strength, in a manner reminiscent of spinodal 
decomposition |3- Note that the presence of free par- 
ticles coexisting with the immobile cluster structure is 
consistent with the presence of charge interactions in the 
poly (methyl methacrylate) system [l2|. 

Fourth, the structural measurements indicate that the 
dramatic slowing down of dynamics as the immobile clus- 
ter phase envelope is traversed could conceivably be ex- 
plained by a combination of two possibilities. On the 
one hand, as the ratio of free to immobilized particles 
decreases (c.f. Fig. Elb)) the mean system dynamics 
becomes increasingly skewed toward localization. On 
the other hand, as pair attractions are increased, par- 
ticles contained in the immobile clusters are themselves 
increasing bound to neighbors and thus localized |27j . 
Analysis of the full measured van Hove self-correlation 
functions (c.f. Fig. [7| indicates that both mechanisms 
are a component of the observed changes in dynamics 
(data not shown). The alternative hypothesis that the 
slowing down of dynamics in the immobile cluster phases 
is due to the quenching of cluster diffusion 56] is not 
an explanation of the Fig. [7| observation. That is, the 
number-density fluctuations of Fig. H3 show that clus- 
ter dimensions are larger than the acquired CLSM image 
volumes. Consequently, our particle-linking algorithms, 
which measure motion relative to the image volume cen- 
ter of mass, are not sensitive to coordinated motions of 
this kind. Instead, the Fig. [7| results report a slowing 
down of the internal dynamics of the immobilized clus- 
ters. 

As pair attractions increase, at a critical c/c*, both 
measures of structure (number density fluctuation and 
contact number) undergo a reversal. We identify this 
transition with a true gelation transition. Structure and 
dynamics become increasingly insensitive to changes in 
pair interaction strength in the gel phase. Interestingly, 
a characteristic structural length scale emerges in the 
system. This scale, about 10 particle radii, is appar- 
ent from the length-scale dependence of number density 
fluctuations (c.f Fig. E{b)). The existence of long-range 
structural heterogeneity in colloidal gels has generated 
considerable controversy. By means of scattering, Ra- 
makrishnan et al have found that depletion gels display 
structural heterogeneity while thermoreversible gels do 
not ^lj. O n the other hand, Varadan and Solomon 
quantified long-range structural heterogeneity in a dif- 
ferent thermoreversible gel system ^1J. In the present 
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FIG. 10: (Color online) (a) Normalized difference between van Hove self-correlation function measured and the Gaussian 
estimate for the resulting mean-squared displacement for c/c* = 0.33. xi represents the point where normalized difference 
between observed and predicted trajectories first becomes negative and X2* represents the trajectories that move past this 
barrier such that there are again more trajectories observed than predicted by a Gaussian, (b) xi*, as a function of c/c*. (c) 
X2* as a function of c/c*. 



report, measurements on well-characterized model sys- 
tem are instructive: the results show that, depending 
on the strength of pair attractions, both kinds of be- 
havior can be found in a single system. That is, at 
low pair attraction strength, long-range structural het- 
erogeneity is apparent in the immobilized cluster phase. 
At high pair attractions, structural heterogeneity is lim- 
ited to intermediate scales (~ 10a). On large scales, 
the gel system returns to a magnitude of number den- 
sity fluctuations that is approximately consistent with a 
homogeneous fluid structure. What remains for future 
work is to resolve how the system selects the particular 
length scale (in this case, ~ 10a) for intermediate range 
structural heterogeneity. A fundamental understanding 
of how this length scale develops in gels could be fruit- 
fully applied to materials development for applications 
involving membranes, supports and scaffolds with par- 
ticular network/void spacing. The existence of this ad- 
ditional structural scale also has strong implications for 
both the linear j5?J and non-linear [35| rheological re- 
sponse of colloidal gels, since predicted properties vary 
by many orders of magnitude in response to structural 
heterogeneity. 

Over a broad-range of pair attractive strengths, the 
depletion system displays strong departures from the 
Gaussian form of the van Hove self-correlation function. 
This behavior is consistent with the observed long-time 
plateau of the single-particle mean-squared displacement 
for both the immobilized cluster and gel phases. To con- 



clude the paper, we interrogate the Fig. 1101 distributions 
to address the possibility of dynamical subpopulations in 
the depletion system. In particular, Kob et al |5^| de- 
scribe a procedure to assess dynamical heterogeneity in 
glassy systems that can easily be applied to gel systems 
as well. The method, applied to the depletion gel sys- 
tem, is shown in Fig. HOf a). We see that, relative to 
the Gaussian reference, there is a surplus of trajectories 
of displacement x(t) < xi*(t), where xi*(t) is defined as 
the first x where G *^ x, *\ , G ^- X,t ' > becomes negative. This 
observation is consistent with localization. That is, an 
excess of trajectories are limited to x < xi*(t). Indeed, 
in the gel phase, Xi*(t) is approximately the range of at- 
tractive depletion potential. Localization on this scale is 
predicted by mode coupling theories of gelation [2]1 • 

If localization were the only consequence of the slow- 
ing down of dynamics in depletion gelation, then xi*(t) 
would be the only dynamical length scale. However, Fig. 
ITU1 suggests an additional dynamical length scale x 2 *(t). 
That is, in both the immobilized cluster and gel phases, 
the subpopulation of colloids with displacements x(t) > 
X2*(t) is anomalously large relative to the corresponding 
Gaussian distribution. This observation suggests that an 
additional mechanism besides the localization predicted 
by mode coupling theory underlies the systems dynamics. 

This region at large displacements of positive devia- 
tions from Gaussian behavior typically comprises about 
1-3% of the systems' total trajectories. In glasses, behav- 
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ior of this kind is an indicator of dynamical heterogeneity, 
the non-exponential decay of time correlation functions 
due to the superposition of individual exponential relax- 
ations with different time constants. Dynamical hetero- 
geneity has been extensively studied in attractive systems 
by means of simulation [28j : however, to our knowledge 
Fig. ^| is the first experimental report of its presence in 
attractive systems. Thus, the results represent a useful 
confirmation that the complex population dynamics of 
gels observed in simulations is indeed consistent with the 
behavior of real systems. 

To conclude, we examine hypothetical origins of the 
observed dynamical heterogeneity. First, in the immobi- 
lized cluster phase, one contribution might be the simple 
superposition of the slow internal dynamics of the clus- 
ters with the fast diffusion of free particles co-existing 
with the clusters. However this explanation is not appli- 
cable to the gel phase because it lacks a significant popu- 
lation of free particles. Alternatively, activated hopping 
transitions are a driver of rare, large displacement dy- 
namics in glassy systems and could hypothetically arise 
in the attractive system studied here |29j| . In addition, 
Puertas et al [3(j discuss the possibility of a collective 
structure origin to dynamical heterogeneity. That is, 
the emergence of structural heterogeneity, such as asso- 
ciated with incipient spinodal decomposition or equilib- 
rium clustering 59, 60], could be correlated with dynam- 
ical heterogeneity. With regards to this last hypothesis, 



we note the coincidence of the system's collective struc- 
ture transition (Fig. EJc)) and the transition in the dy- 
namical heterogeneity length scale X2* (Fig. llOf c)). In 
addition, the strength of pair attractions appears to af- 
fect collective structure in the immobile cluster and gel 
phases in a way that is similar to its effect on dynamical 
heterogeneity. That is, collective structure and dynami- 
cal heterogeneity are much more sensitive to the strength 
of pair attractions in the immobilized cluster phase than 
in the gel phase. Thus, although we cannot at this time 
definitively resolve among the listed hypotheses, we note 
that the comparison of the features of Fig. |SJc) and Fig. 
HOl c) are highly suggestive of a role for collective struc- 
ture in mediating dynamical heterogeneity in attractive 
gel systems at intermediate volume fractions. 
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